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Abstract 

We consider supersymmetric unification based on the gauge symmetry 
SU(3)c x SU(3)l x SU(3)r supplemented by an anomalous U{1)a gauge 
symmetry. Realistic fermion masses and mixings are realized, including 
large mixings in the neutrino sector. We also consider the supersymmetric 
flavor problem, gauge coupling unification and proton decay. The dominant 
proton decay mode is expected to hep — > e + 7r° and the lifetime is estimated 
to be ~ 10 34 - 10 35 years. 
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1 Introduction 



Supersymmetric grand unified theories offer a particularly elegant scenario for 
unifying the strong and electroweak interactions. Recently, realistic models based 
on 50(10) 0, |3| and E 6 ^, |J have been proposed in which an anomalous U{1)a 
gauge symmetry]^] plays a critical role. The anomaly is cancelled via the Green 
Schwarz mecanism|0], and the resulting phenomenology has several attractive 
features. In particular, all interactions allowed by the symmetries are included 
in the discussion (there can be undetermined order unity coefficients accompany- 
ing the interactions. These models naturally resolve the doublet- triplet problem 



using the mechanism of reference ||10|| (see also [111], 12, [13|]). Realistic pattern 



of qaurk and lepton mass matrices, including large neutrino mixings [14]], are real- 
ized, using the Froggatt- Nielsen (FN) mechanism [|15|. The anomalous U(1)a also 
helps explain the hierarchical symmetry breaking scales and the masses acquired 
by the superheavy particles, and leads precisely to the minimal supersymmetric 
standard model (MSSM) at low energies. Even though the gauge couplings unify 
in these schemes slightly below the usual GUT scale 2 x 10 16 GeV, dimension 

five proton decay is sufficiently suppressed , and the decay p > e + n via 

gauge mediated dimension six operators may be seen in the near future. Finally, 
in these models the cutoff scale is lower than the Planck scale Mpi anc ^ and the /z 
problem is also resolved. 

However, these models require two adjoint Higgs fields to realize DT split- 
ting, which is not so easily realized in the framework of superstring models. In 
this paper we examine the application of the above approach to grand unified 
theories with semi-simple unification whose symmetry breaking to MSSM (the 
minimal supersymmetric standard model) does not require adjoint Higgs field. 
A particularly attractive example is provided by the gauge symmetry SU(3) 3 = 
SU(3)c x SU(3)l x SU(3)r which is a maximal subgroup of E 6 , and which 
arises as an effective four dimensional symmetry from the compactification of the 
E$ x E 8 heterotic superstring theory on a Calabi-Yau manifold ||16||. Phenomenol- 
ogy based on SU(3) 3 has been extensively discussed in the past [IT], [TS|, [B|. Our 
goal here is to apply the techniques of |l|-[f| to gauge symmetry SU(3) 3 and elu- 
cidate the most important consequences. In particular, we will show how realistic 
fermion masses and mixings are obtained, including large mixings in the neutrino 
sector (with the exception of U e3 ). We also consider SUSY breaking and flavor 
changing neutral currents, gauge coupling unification and proton decay. While 
dimension five proton decay turns out to be heavily suppressed, dimension six 
operators yield a proton lifetime of ~ 10 34 — 10 35 years for the decay channel 
p — > e + 7i°. 
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2 Matter sector 



The matter sector has essentially the same structure as the Eq model 0, with 
the 27 of E 6 given in terms of SU(3) C X SU(3) L x SU(3) R as 

27 -(3, 3,1) + (3, 1,3) + (1,3, 3). (2.1) 

Three 27-plets (i = 1, 2, 3) are introduced, and the Yukawa interations contain 
appropriate powers of the VEV of FN field (©) = AA, which has an anomalous 
U(1)a charge 9 = —1, namely 

X^+^+^i^fj®. (2.2) 

Here $ is a Higgs field, and throughout this paper we use units in which the 
cutoff A = 1, and denote all superfields by uppercase letters and their anomalous 
U(1)a charges by the corresponding lowercase letters. Using the definitions of the 
fields Q(3, 2) i , f/ c (3, 1)_ § , D c (3, 1) i , L(l, 2)_i , £ c (l, 1) 1; N c (l, 1) , L c '(l, 2) i , 
L'(l,2)_i, D'(3,l)_i, ^(3,1)1, 5(1,1), and their conjugate fields under the 
standard model (SM) gauge symmetry, the fields (3, 3, 1), (3, 1, 3), and (1, 3, 3) 
under SU(3) C x SU(3) L x SU(3) R are 

(3,3,1) -> Q + ZT (2.3) 
(3,1,3) -> C/ C + D c + D c/ (2.4) 
(1,3,3) -> L + E c + N c + L ' + L' + S. (2.5) 

For future reference, under the breaking E e — > 50(10), 

27 -» [g + [7 C + E c + £> c + L + iV^ + jL^ + L 7 + D' + D c/ ] j + ^. (2.6) 

16 10 1 

Since D'(L C/ ) can acquire superheavy mass by combining with a linear combina- 
tion of D c and D c ' (L and II) after the breaking SU(3) C x 5f/(3) L x Sf/(3) fl 
to the SM gauge group, the remaining massless fields form the three generation 
matter content of MSSM. Since the Yukawa couplings are determined mainly by 
the anomalous U(1)a charges of the massless fields, we would like to know which 
of the fields among D\ and Df (Li and L[) (i = 1, 2, 3) are massless. 

To this end, we discuss how SU(3) C x SU(3) L x SU{3) R breaks to MSSM. 
We introduce the Higgs fields with non-vanishing VEVs $(1,3,3), $(1,3,3), 
C(l,3,3), and (7(1, 3, 3). The VEVs | ($) | = | ($)| ~ _\-(0+0)/2 ^reak SU (3) c x 
SU(3) L x 5^(3)^ to SU(3) C x S'C/(2) L x 5^(2)^ x U(1) B -l, while | (C) | = 
I (C) I ~ A-( c+c ~)/ 2 break 5C/(3) C x SC/(2) L x 517(2)^ x C/(1) B _ L to the SM gauge 
group. Here the VEVs are determined by the anomalous U(1)a charges, and the 
reason is roughly as follows (We will explain how to determine the VEVs later): 
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1. Since the interactions are determined by the anomalous U(1)a charges, the 
VEV of the gauge invariant operator O with negative charge o is determined 
as (O) ~ A"°. 



2. For O = the VEV becomes ~ A" 



3. Since the D-flatness condition requires | | = | ($) |, we obtain | ^$ 

The massless modes can be determined from the superpotential 

W\ = A^^ + **i*j$ + A^+^+^^C, (2.7) 

and the VEVs ($) ~ A^W'+'W and (C) ~ A~^ (c+5) . Here, for simplicity, we have 
assumed the Eq like charge assignment in the matter sector, but in principle, we 
can assign these charges without respecting E§ symmetry. The mass matrices of 
D c ^ and D' {L® and L cl ) are obtained from 

h h h I\ 1*2 I3 

^yl)\+ip3+r ^ip2+*l>3+r y2^3+r ^1+^3 \i>2+i>3 \^3 J 

(2.8) 

(c-c-0 + 0), (2.9) 



I' 3 



where I = L, D c and we have defined the parameter r as 

_ 1 
r = -I 

which we use frequently in the following discussion. Note that the mass matrices 
are determined by the anomalous U(1)a charges. Therefore the massless modes 
are also determined by the charges. As discussed in Ref . , as long as we neglect 
the cases with vanishing coefficients arising from some SUSY based mechanism, 
the main components of the massless modes can be obtained as follows: 

1. Vi ~i>3 < r : (h,I 2 ,h)- 

2. 0<r <^i-^3 : {h,I[,h). 

3. V>3-V>l<r<0: {h,I[J' 2 ). 

4. r<^-^:[I[,I' 2 ,I' z ). 

The case (J 1; /{, I 2 ) is interesting, because bi-large neutrino mixing angles can be 
realized without tan (3 too small, if we take account of the mixing of subcompo- 
nents. Indeed, the massless modes (1°, 1° , I3) are given by 

J° = Il + x^-^I 3 + \^-^+ r I' 2 + \^-^+ r I' 3} (2.10) 
J 2 ° = I[ + \i>i-i>3-rj 3 + yh-ihj'^ + yPi-^sj^ (2.11) 
J 3 ° = J 2 + \^-i>3j 3 _|_ j/ _|_ ^ 2 -^3+rj^ (2.12) 
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where the first terms on the right-hand sides are the main components of these 
massless modes, and the other terms represent mixing with the other states, J 3 , 
I' 2 and 

The mass matrices for quarks and leptons are obtained from the superpoten- 

tial 

W 2 = X^+^+H^jH, (2.13) 

where #(1,3,3) contains the MSSM Higgs doublets. If we adopt the charges 
■0! = 3 + n, — 2 + n, ifti — n, and h = — 2n, the mass matrices are given by 



Mu = 



Qi 
Q2 



U{ ui ui 
( A 6 A 5 A 3 \ 
A 5 A 4 A 2 



(L C '(H)) , 



m d (m|t 7 - 1 ) 



Q3 \ A 3 A 2 I J 

Df(Ll) D?(L°) Df(Ll) 
Qi{El) I A 6 A 6 ~ r A 5 
Q2(E C 2 ) 



A 5 
A 3 



A 5-r 
A 3-r 



A 4 
A 2 



(2.14) 



(£'(#)>. (2.15) 



Here i] ~ 2 — 3 is the renormalization group factor. Then we can obtain the CKM 
matrix as 

/ 1 A A 3 \ 



U, 



CKM 



(2.16) 



A 1 A 2 
V A 3 A 2 1 J 

which reproduces the experimental results if we take A ~ 0.2. Since the ratio 
of the Yukawa couplings of the top and bottom quarks is A 2 , a small value of 
tan/3 ~ mt/mb ■ A 2 is predicted by these mass matrices. 

The Dirac neutrino mass matrix is given by the 3x6 matrix 





Si S2 


S 3 




N c 2 


iV 3 c 


LI 




A r+3 


A 6 


A 5 


A 3 


L°2 


A 6 A 5 


A 3 


A 6-r 


A 5-r 


A 3-r 


LI 


{\ r+5 A r+4 


y+2 


A 5 


A 4 


A 2 


simply express as 














( 


A 4 


A 3 


A 



M N = ( \ r+2 A 2 ) 



A 4-r 

V A 3 



A 3-r 

A 2 



(L C, (H)) V , 



A l-r 
1 



(2.17) 



(L«(H))r,. (2.18) 



The 6x6 matrix for the right-handed neutrinos (Si,i = 1,2,3, and N%,k = 1,2,3 
) is obtained as 



M, 



R 



XPi+^+^SiSi M + A^ + ^ +5+ ^^iV fc c ($\ (c) , 
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I \4>~4> \{i-4>+c-c)/2 \ / A 6 A 5 A 3 \ 

- ( x~ a- ) ° ( j: j; a ; J (22o > 

from the interactions 

A^ + ^ +2 ^* i ^$$ + A^+^+^^C' + A^ i+ ^ +2 %# i C'C', (2.21) 

by developing the VEVs ($) ~ and ~ A~5( c + C ~). Using the seesaw 

mechanism p0| , we obtain the neutrino mass matrix 

/ A 2 A 2 " r A \ 
M v = MnM^MJi = A 4_2n+C_e A 2 - r A 2 - 2r A 1_r (L C, (#))V- (2.22) 

V A Al ~ r 1 / 

Then, we finally obtain the Maki-Nakagawa-Sakata matrix as 

/ 1 A r A \ 
U U m = \ X r 1 A 1 ^ . (2.23) 




If we adopt r = 1/2, namely, 

c -c = 0-0 + l, (2.24) 

bi-large neutrino mixing angles are realized because of A 1 / 2 ~ 0.5 as within the 
SO (10) model in Ref. |J and E% model in Ref. [|J. Moreover, it predicts 
U e3 ~ A. Future experiments can see whether there is such a large U e3 just below 
the CHOOZ upper limit U e ^ < 0.15 [ 2TQ or not. For the neutrino masses, the 
model predicts vn Vti jm VT ~ A, which is consistent with the most probable LMA 
MSW solution for the solar neutrino puzzle. [14| 
If we define 

Z = 0-0 + 2n-lO, (2.25) 
the neutrino mass matrix is given by 

/ A 2 A 15 A \ 
M v = A~ (5-M) A 1 ' 5 A A ' 5 (L C '(H)) 2 t] 2 , (2.26) 
V A A 05 1 J 

where we have used the relation ( 2.24 ). The paramter I can be determined by 

A . = ^™v, (2 . 27) 
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We are supposing that the cutoff scale A is in the range 5 x 10 15 GeV < A < 
10 20 GeV, which allows —3 < I < 2. If we choose I = —2, the neutrino masses are 
given by m Vr ~ A" 3 (L C '(H)} 2 r/ 2 /A ~ m v J\ ~ m^ e /A 2 . If we take r? (L c, (#)> = 
100 GeV, A ~ 2 x 10 16 GeV and A = 0.2, then we get m Vr ~ 6 x 10~ 2 eV, 
m v ~ 1 x 10 -2 eV and m Ue ~ 2 x 10 -3 eV. These values are pretty much 
consistent with the experimental data for atmospheric neutrinos and with the 



large mixing angle (LMA) MSW solution for the solar neutrino problem. 22 



3 SUSY breaking and FCNC 

Let us now discuss SUSY breaking. Since the anomalous U(1)a charges depend on 
flavor to produce the hierarchy of Yukawa couplings, generically non-degenerate 
scalar fermion masses are induced through the anomalous U(1)a -D-term. Under 
the Z?6-hke charge assignment of the matter sector, the SUSY contribution to 
K° — K° mixing is naturally suppressed as in Eq GUT scenario. The essential 
point is that the anomalous U(1)a charge of Df becomes the same as that of 
D?, , because the field ~ D\ and Df ~ Df arise from a single field of Eq. 
Since the constraints from K° — K° mixing to the ratio 5 = A/m 2 , where A is 
the mixing of sfermion mass matrices and rh is the average of sfermion mass, are 

\M8&)ll(6? 2 )rr\ < 2.8 x 10- 3 (^M^) ' (31) 

\Re(5°) LL \, \Re(6°) RR \ < 4.0 x 10" 2 (^^) , (3.2) 

and the former constraint is much stronger than the latter constraint, suppression 
of (8® 2 )rr makes the constraint on the SUSY breaking sector much weaker. Here 
we use the notation (Sfj)xY, where F = U, D, N, E, the chirality index is X, Y = 
L, R, and the generation index is i,j = 1, 2, 3, as defined in Ref. p5 |. As in the 
usual anomalous U{1)a scenario, A can be estimated as 

(^xx~\ ]fi - f *K\fi-fj\)(DA). (3.3) 

Therefore, the coincidence of the anomalous U(1)a charges of Df* and D^ leads 
to the suppression of (A^ 2 )rr- This weakens the constraints from K°—K° mixing. 

To see how weak the constraints become in our scenario, we fix the SUSY 
breaking sector as follows. At the cutoff scale, we adopt the common gaugino 
mass Mi/2 and D-term of anomalous U(1)a gauge symmetry Da as no n- vanishing 
SUSY breaking parameters. Then, the scalar fermion mass squared at low energy 
scales is estimated as 

m%~ fiRMl /2 + VF Ml /2 , (3.4) 
where r]p is a renormalization group factor and 

R S ™. (3.5) 

1VI l/2 
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The constraint ( |3.2| ) for (5^)ll is rewritten 



M 1/2 > 1.25 x 10 % ^, ^„\ 3/9 (GeV). (3.6) 



Though the main contribution to (S^rr vanishes, through the mixing in Eqs. 
(|2.10|) and Q2.ll]) , (5® 2 )rr is estimated as 

{o 12 )rr ~ A2 . (3.7) 

vd r +ipiR 



From Eq. ( [3.1|) for y (5i2)ll(<5i2),R-R> the constraint on the gaugino mass Mi/ 2 is 
given by 

M in > 1.8 x iff . (3.8) 

On the other hand, the \i — * e7 process gives 

l(*f 2 M, K^M < 3-8 x 10~ 3 (^p^ , (3-9) 

where m; is the average mass of the scalar leptons. This constraint is rewritten 

.( A(^-^)fl) 1 / 2 



^ > 1-6 >< 10" ^:^ (GeV). (3.10) 



Taking the values = 9/2, ^ 2 = 7/2, t] Dl ~ 77^ ~ 6 and tje r ~ 0.15, the 
rough lower limits on the gaugino mass are in Table I. 



Table I. Lower bound on gaugino mass Mx/2 at the GUT scale (in GeV). 



R 


0.1 


0.3 


0.5 


1 


2 




17 


43 


61 


87 


105 


\J {&vi)ll{$x2)rr 


78 


191 


268 


373 


437 


m)RR\ 


431 


304 


221 


161 


116 



Note that in some range of R, the fi — > c~) process gives the severest con- 
straint among the FCNC processes. |2B] Therefore, the lepton flavour violating 
processes p6|, |27| might be seen in the near future. 



4 Higgs sector 

In addition to the Higgs with non-vanishing VEVs $, $, C, and C, we introduce 
C'(l, 3, 3), C"(l, 3, 3) with vanishing VEVs and several singlets S and Z in order 
to give superheavy masses to these Higgs fields. The Higgs content is 
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Table II. The typical values of anomalous U(1)a charges are listed. 





non-vanishing VEV 


vanishinj 


y VEV 


(1,3,3) 


$(0 = -3) C(c= -3) 


C"(c' : 


= 0) 


(1,3,3) 


$(0 = 2) C(c= 1) 


C"(C' : 


= 4) 


1 


Z i (z i = -l)(i = 1,2,3) 






(3,3,1) 




Qh{qi 


= 1) 


(3,3,1) 




Qh(qi 


= 0) 


(3,1,3) 




Qn{qr 


= 1) 


(3,1,3) 




Qn(qr 


= 0) 



Here the Higgs field H is contained in <3> as in the E 6 case. The Higgs fields Ql, 
Ql, Qr, and Qr are introduced only for realizing the same Kac-Moody levels of 
the three SU (3) gauge groups and they do not play any other role in the following 
argument. 

In this model, the singlet composite operator plays the same role as the 
FN field O. The D-flatness condition for the anomalous U(1)a gauge symmetry 
is 

D A = gA (e+m 2 +m 2 ) = o, (4.1) 

where £ 2 is the parameter of the Fayet-Illiopoulos D-teim. Since the D-flatness 
conditions of SU(3) L and SU(3) R require | ($) | = | |, the D-flatness condi- 
tion for the anomalous U{1) A gauge symmetry is rewritten 

D A = g A (e + (<P + ^m 2 ) =0. (4.2) 

Thus we obtain £ 2 + (0 + 0)|$| 2 = 0, namely, | ($) | = = £. In this case, 

since plays the same role as O, the unit of hierarchy becomes = A ~ 

which is different from the usual case in which the FN field is just a singlet field 
and (6) = £. It means that even if £ has a milder hierarchy, the unit of 
hierarchy becomes stronger. Using gauge rotation and D-flatness condition for 
SU(3)l x SU(3)r gauge symmetry, the VEV can be taken as 

/0 \ 

| ($) | = | ($\ | = , (4.3) 

\o o x-hio+ft) 

which breaks SU(3) C x SU(3) L x SU(3) R into SU(3) C x SU(2) L x SU(2) R x 
U{1)b-l- In order to determine the VEVs of the other Higgs fields, we examine 
the following superpotential 

W = W C " + W e > + W NV , (4.4) 

where Wx denotes the terms linear in the field X, which has vanishing VEV, and 
Wnv includes only the fields with non-vanishing VEVs. From the superpotential 

W NV = $ 3 + ¥C, (4.5) 
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(l\C)) = (L C '(C)) = is obtained. The vacuum is ((S(C)) = 0, (N C (C)) ^ 0) 

or ((S(C)) ^ 0,(n c (C)) = 0) if (Cc) ^ 0. We are interested in the first 

vacuum {(S{C)) = 0, (N C {C)) + 0). 

The superpotential W c > and Wq, are given by 

W c > = A C ' + ^C"<5(1 + X E+C CC + \ Z ^Z % Z 3 + X Z 'Z Z ) 

+\ C ' +5 CC(1 + X Zl Z t ) (4.6) 
W & = A c "' + *C"<l>(l + A Zl Z i ) + A c "' +c C"C , (l + A 21 ^). (4.7) 

Here we neglect ($<3>) 2 for simplicity, but the effect is critical. After developing the 
VEVs, the above interactions do not respect SU(3)l x SU(3)r gauge symmetry. 
For example, the coefficient of N C (C)N C (C) is different from that of L{C')L(C). 
This is important to align the VEVs and to give superhevay masses to these fields. 
The F-flatness conditions F S (c') — ^s(C") = -^v c (C") — Fn c (C) = determine four 
VEVs (Cc) ~ X^ c+5 \ and (Zi) ~ \~ z *(i = 1,2,3). Then all the VEVs are 
determined by the anomalous U(1)a charges. 

We now examine the mass spectrum of the Higgs sector. The mass matrix 
Ml for L and L is obtained from the interactions 







TO 
Ll§ 


TCI 

L c 


L% 


L'c 


TCI 

L,Q, 


Tl 

L,Q, 




Lc 






T> 


( 






























L'c 



















cc 























$3 













C"C$ 




TO 

L c 










3> 2 C 


§c 2 


cc 




c 2 § 


c 


C"C 2 




Tl 

Li c , 










C"$ 


cc 





cc 











TCl 
LiQ, 




C"$ 


cc 


C"$ 2 


cc$ 


cc 


C" 2 $ 


c"c$ 


cc 2 


C 2 C 2 























CMC 








C"$ 




L c 



















c%c 2 








C"C 




Lc 


\ 








<§ 2 C"C 


Wcc 





c%cc 


C"$ 


cc 


cc 


) 



(4-8) 

It is obvious that the linear combination of L' $ and L' c , and that of and L' c are 
massless and they form the doublet Higgs fields of MSSM . Lq> and Lq, are eaten 
by the Higgs mechanism in breaking SU(3)l x SU(3)r into SU{2) L x SU(2) R x 
f(l)s-L- The mass spectrum of the remaining fields becomes A c+C ', A c+C ', A c+C ', 
A c '+ 5 , A c ' +C ", A c '+ 5 , and A 2< * + ^-<«. 

The mass matrix Me for i? c and E c is obtained from the interactions 





H 


E°- 
^c 


% 




( 







E c c 











E c , 


{c$ 





cc) 
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The fields E c c and are eaten by the Higgs mechanism in breaking SU(2) R x 

U(1)_b-l into U(1) Y - The mass spectrum of the remaining fields is A c ' +<?i and 

The mass matrix for the fields D c , D cl , D c , D c ' is obtained from the in- 
teractions 



D% D% D' n 

Qr Qr Ql 

D c Qr ( QrQr Qr^QrC \ 



U Qr 

D 'o 

Ql 



(4.10) 



QrQr 

\QlQrC Q l Qr§ Q l Q l J 

The mass spectrum becomes A 5r+9r , A 5r+9r , and \^ +qi . 

The mass of the fields Q and Q is obtained from the interaction QlQl as 
\qi+h_ The mass of the fields U c and U c is obtained from the interaction QrQr 

as 

By the above argument, the mass spectrum of superheavy particles are deter- 
mined only by the anomalous U(1)a charges, so we can examine whether coupling 
unification is realized or not. Before going to the discussion in the next subsec- 
tion, we define the reduced mass matrices Mj by getting rid of the massless modes 
from the original mass matrices Mj. The rank of the reduced matrices in our 
semi-simple model are ?q = f Uc = 1, f E c — 2, f L — 7 and f D c = 3. It is useful to 
define the effective anomalous U(1) A charges: 

Xj = i + I A 0, x T = J-^,(x = l',l c ',d c ',d'), (4.11) 

1 1 

xi = i + Ac- -A0, xj = i - Ac+ -A0, (x = l,d°), (4.12) 

xi = i, xj = i, (x = q, u c , e c ), (4-13) 



2(C-C) 



where I = §,C,C',Q L ,Q R (i = (p,c,c',q h q r ), A0 = |(0-0), and Ac 
The determinants of the reduced mass matrices are estimated by simple sums of 
the effective anomalous U{1)a charges of massive modes: 

det M Q = X^L +q QL = \qi+qi (4.14) 
detMc/c = X^Qr+iQr = \^+ir ( 415 ) 

detM £c = A e l +e * +e "c' +e c' = A^ +5 ' +c ' (4.16) 

detM^c = \ d, QL +d 'Q L + %R +d QR +d QR +d QR = A ffl+«+2(5r+tt0 (4.i7) 

det M L = A z c+^+^+'S'+^'+ r c+^+ z c+%'+'"g+'c'+ r e' + ^+«c' 

_ ^3(c+c+c'+c')+2<£-A</>_ (4.18) 

Then all the elements of mass matrices are estimated by simple sum of the effec- 
tive charges of superheavy particles if they are not vanishing, and the determi- 
nants of mass matrices are also determined by simple sum of the effective charges. 
We will use this result in calculating the running gauge couplings. 
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5 Coupling unification 



In this section, we apply the general discussion on the gauge coupling unification 
in Ref. to our scenario. The pattern of the breaking of the gauge symmetry 
in our model is as follows. At the scale A$ ~ A - ^"^" 2 , SU(3) 3 is broken into 
SU(3) C x SU(2) L x SU(2) R x U(1) B _ L . The SU(2) R x U(1) B -l is broken into 
U(1) Y at the scale A c ~ A~ (c+5)/2 . We base our analysis on one loop renormaliza- 
tion group equations. The conditions of the gauge coupling unification are given 
by 

a 3 (A) = a 2 (A) = -Q!y(A) = ai(A), (5.1) 

where a^ 1 (A$ > /i > Ac) = (A$ > fi > Ac) + |a^i L (A$ > [i > 
Ac), cti 1 ^ > A$) = ia-^/i > A$) + ha^Kfi > A$), and a^/i > A$) = 

2 

oi^ L {ii > A$) + > Ac). Here ax = |^ and the parameters gx(X = 

3,3L,3R,2, R, B — L,Y) are the gauge couplings of SU(3)c, SU(3)l, SU(3) r , 
SU(2) L , SU(2) R , U(1) B -l and U(1) Y , respectively. 

Using the fact that the three gauge couplings of the minimal SUSY standard 
model (MSSM) meet at the scale Ac ~ 2 x 10 16 GeV, the above conditions for 
gauge coupling unification are rewritten 

>Hi>^H^y>(t)-Hi) <"> 

= ^ ln (^) + E ' Ai2 ' ln (dS;)- 21n (^) (5 - 3) 

= 6 3 ln(^) + E,A63,ln(^), (5.4) 

where (&i, 62,^3) = (33/5,1,-3) are the renormalization group coefficients for 
MSSM and Ab a j(a = 1,2,3) are the corrections to the coefficients from the 
massive fields I = Q + Q,U C + U c , E c + E c , D c + D c , and L + L. The second 
last term in Eq. (|5.2| ) is from the breaking SU(2) R x U{\)b-l — > U(l)y by 
the VEV (C), and the last terms in Eqs. ( |5.2| ), ( |5.3| ) are from the breaking 
SU(3) L x SU(3) R -> SU(2) L x SU(2) R x U(1) B -l by the VEV ($). Since all 
the mass matrices and the symmetry breaking scales appearing in the above 
conditions are determined by the anomalous U(1)a charges, these conditions can 
be translated to the constraint on the effective charge and the cutoff scale, 

AgA^ (5< ^ +6(c+c " )} , (5.5) 
A G Ai (7< ^ +6(c+5+c ' +5 ' )} , (5.6) 

AG ,^gL(25^+19^+30( C +e)+18(c'+5')_ ( 5 _ 7 ) 

A naive calculation leads to the relation between the charges, 

590 + 410 + 42(c + c) + 54(c + c') ~ 0, (5.8) 



ai(A) 


= MA) - 


-> A 


a 2 (A) 


= 013(A) - 


-> A 


ai(A) 


= a 3 (A) - 


-> A 
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which is difficult to satisfy in our scenario. However, careful calculation shows 
that gauge coupling unification is possible, though somewhat larger ambiguities 
of order one coefficients are required than in a simple group unification. Actually 
with the typical charge assignment in Table I, the coupling unification is realized 
as in Fig. 1, using the ambiguities of order one coefficients A < y < A . 




2.5 5 7.5 10 12.5 15 log/i(GeV) 

Figure 1: Here we adopt A = 0.22, a^{M z ) = 59.47, a^{M z ) = 29.81, 
aj^Mx) = 8.40, the SUSY breaking scale mss ~ 1 TeV and the anomalous 
U(1) A charges: (^1,^2,^3) = (9/2,7/2,3/2), = -3, = 2, c = -3, c = 1, 
z = —1, d = and c' = 4. Using the ambiguities of coefficients A < y < A -1 , 
three gauge couplings meet at around 10 16 GeV. 



Note that the unified gauge coupling can remain in the perturbative region in 
SU(3) 3 model, which is different from the Eq model. This is because the Higgs 
sector in SU (3) 3 unification is much simper than that in E 6 unification, since we 
do not have to introduce adjoint Higgses to realize the doublet-triplet splitting. 

The cutoff scale tends to be lower than the Planck scale. Indeed, the cutoff 
is taken as 10 16 GeV in Fig. 1. Since the cutoff scale is so low, we have to take 
care of proton decay via dimension five operators [3(|, which are obtained from 

A^+^+^+^+^^fctyf (5.9) 

by developing the VEV ~ A^^+'W. The coefficients are suppressed not only 
by the usual small Yukawa factor but also by the suppression factor \ in+ \i't > -'t > = 
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A 8 ' 5 . Even if we take the cutoff A ~ 10 16 GeV, the 'effective' colored Higgs 
mass is around A~ 85 A ~ 10 22 GeV, which is much larger than the experimental 
bound of 10 18 GeV. Thus, proton decay via dimension five operator is adequately 
suppressed. 

On the other hand, proton decay p — > e + ir° via dimension six operators from 
the Kahler potential 

if^iflM^i, (5.10) 

which are allowed by the symmetry in our scenario by taking the unification scale 
Au as the cutoff A, may be seen in future experiments. If we roughly estimate 
the lifetime of proton using the formula in Ref. JHJ and the recent result of the 
lattice calculation for the hadron matrix element parameter q ||33| , we find 



r Ap e V) ~ 4.5 x 10" (Mj year, ,5,1) 

This estimate, albeit a rough one, provides a strong motivation for continuing 
the proton decay search. 




6 Discussion and Summary 

Besides SU(3) 3 , E e has the other maximal semi-simple subgroups SU(6) x SU (2) l 
and SU(6) x SU(2)r^%. The matter sector can be applied to these subgroups 
in a straightforward way. However, in the Higgs sector, it is difficult to realize 
the situation in which only one pair of doublet Higgs is massless. It is difficult to 
make the partner of the doublet Higgs massive, while keeping the latter massless. 
On the contrary, in SU(3) 3 gauge symmetry, since the partners of the doublet 
Higgs L, E c and A^ c are absorbed by the Higgs mechanism, one pair of doublet 
Higgs can be massless. 

In the typical charge assignment, the charges of the matter sector respect the 
E 6 symmetry, while those of the Higgs sector do not. It is difficult to respect E 6 
symmetry in the Higgs sector without additional massless fields other than the 
fields in the MSSM. 

By introducing singlet fields, we can build models with integer Kac-Moody 
level. For example, in addition to the fields in Table I, we introduce a singlet 
field with charge 10, one with charge -8, 43 singlet fields with charge 3/2, and 62 
singlet fields with charge 1/2. Then using the relation [] 



1 C a = Trc a T(R) Qa- Here T(R) is the Dynkin index of the representation R, and we use 
the convention in which T (fundamental rep.) = 1/2. 
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where k A and k a are Kac-Moody levels of U(1)a and SU(3) a (a = C, L, R), these 
Kac-Moody levels can be calculated as 

k A = 4, k c = k L = k R = 2. (6.2) 

Note that introducing the singlets with charge 10 and -8, the \x problem is solved 
by the mechanism proposed in Ref. [[|. 

In our model the difference between the mass matrices of down-type quark and 
the charged leptons are realized because the matrices are from different Yukawa 
interactions. However, if this model is regarded as the low energy theory of 
E 8 x E s heterotic superstring theory, we have to break the gauge symmetry E e into 
SU (3) 3 . Since the matter sector respects E 6 symmetry, it is natural to expect that 
the Yukawa interactions also respect it. In order to realize the different Yukawa 
interactions, we have to implement the breaking. In the brane world scenario, 
there is an interesting mechanism to break the gauge symmetry 0. However, it 
seems difficult to realize the breaking in the Yukawa coupling of matter which 
resides on the brane. To enforce the E$ breaking some of the matter must be in 
the bulk, where the E G gauge symmetry is not respected. 

In this paper, we have proposed a realistic semi-simple unified theory with 
SU(3) 3 gauge group. Since generic interactions have been introduced, we can 
define the model by the anomalous U(1)a charges. Large neutrino mixing an- 
gles can be realized in the model. Moreover, the FCNC process is automatically 
suppressed. The half integer charges of the matter sector automatically play the 
role of .R-parity. The model has the same matter structure as in Eq model |3J], 
but different and simpler Higgs structure. Actually we do not need the adjoint 
Higgs fields to realize doublet-triplet splitting. And the gauge couplings at the 
cutoff scale can be in the perturbative region. Note that in the SU(3) 3 model, in 
contrast to SU(5), 5*0(10) or E 6 , the lightest magnetic monopole carries three 
(instead of one) quanta of Dirac magnetic charge^]. This is readily seen by 



noting that one is allowed, in principle, to include non- bifundamental vector- 
like representations such as (1,3,1) + (1,3,1) that, despite their color singlet 
nature, carry fractional(e/3) electric charge. The Dirac quantization then re- 
quires that the corresponding magnetic charge has three units. The number 
density of primordial SU(3) 3 monopoles depends, of course, on the underlying 
cosmological scenario, and should not exceed the nominal Parker bound of about 
10~ 16 cm -2 s -1 sr -1 . The discovery of magnetic monopoles would be a truly re- 
markable event, and measurement of their magnetic charge would allow us to 
distinguish between a variety of unified gauge theories. 

It would be interesting to extend the approach presented here to other semi- 
simple unification schemes. For instance, the gauge symmetry SU (3) 3 with three 
27's of E e can be embedded, in principle, in SU(A) x SU(3) x SU(3)$H\ which 
could be worth pursuing. 
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